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J^ffiTHOD OF PROCESSING A SIGNAL USING AN APPROXIMATE MAP 
ALGORITHM AND CORRESPONDING APPLX^ATT^N^ , 



The domain of the invention is processing of signals 
transmitted in a channel that is or that might be 
disturbed. 

More precisely, the invention relates to a signal 
5 processing method making use of a MAP (Maximum A 
Posteriori) type algorithm. The MAP algorithm is 
presented particularly clearly in the article entitled 
"Statistical Inference for Probabilistic Functions of 
Finite State Markov Chains" by L. E. Baum and T. Pertie, 
10 Ann. Math. Stat. 37: 1554 - 1563, 1966. 

The MAP algorithm, also called the Forward-Backward 
(FB) algorithm, is drawing a great deal of interest at 
the moment. In particular, this algorithm is capable of 
decoding a convolutive code by associating reliability 
15 information with decoded bits, and it can beneficially be 
used in many applications such as applications using a 
"turbo-code" . 
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The following briefly describes the principle of the 
MAP algorithm. For the purposes of this description and 
throughout the remainder of the document, the term "MAP 
algorithm" includes MAP type algorithms, and also 
5 algorithms derived from MAP type algorithms, and 
particularly Log-MAP and Max-Log-MAP type algorithms . 

Diagrammat ically , the MAP algorithm may be 
considered as merging of information (node metrics) 
generated by two modified Viterbi decoders working on 

10 data in the backwards direction. Remember that the 
Viterbi algorithm, frequently used for decoding of 
received signals, consists of searching in a lattice 
representing the coding for the path corresponding to the 
most probable sequence of states, in other words 

15 searching for the path which is the minimum distance from 
the required sequence. This algorithm is described in the 
article by G. D. Forney, "The Viterbi Algorithm", 
Proceedings of the IEEE, vol. 61, No. 3, March 1973. 

The MAP algorithm can be used to determine the 

20 likelihood ratio of the a posteriori probability of 
states at instant k, 0 ^ k ^ K of a hidden Markov string 
with finite states, starting from observations 

= ^ Yo,...^ Yfc}- . A Markov chain is usually defined as a 
collection of consecutive random states chained one 
25 behind the other, in which the changeover from one state 
to another is associated with a probability called the 
transition probability. 

A Markov chain ^X^}- with values in a finite space 
E = can thus be considered, defined by a states 
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transition matrix f] = i^i.j) r with dimension N x N and an 
initial law v = (vi) defined by Vi = Projb[Xo = i], Vi e E. 
According to the terminology conventionally used, a 
transition between two states Xk-i and Xk is defined by a 
5 pair of states and is denoted = (Xk-i, Xk) . The matrix Yl 
is a transition matrix between states associated with 
for which the elements nxj are defined by: 
jtij = Projb[Xk = j/ Xk-i = i ]/ f or all i,j G E (1) 
This document considers the case of hidden Markov 
10 models for which noisy observations are available \ Yk}- 
with value in 91*^, collected through a channel without 
memory. In other words, observations ] Yr} are mutually 
independent depending on transitions between states ^a:^ 
and each observation Y;: is only dependent on the 
15 transition ^a: at the same instant k. Consequently, the 
observed signal model can be written in the following 
form: 

Yjc = H(^j,) -f Bj, (2) 
where H{.) represents the transfer function of the 
20 channel without memory defined in E and with value in 91*^, 
^/c = (X/c-i, X;c) is a pair of states at instant k with value 
in E X E, Yk represents an observation at instant k with 
value in di^ and where is added white noise. 

In the context of digital communications, the model 
25 defined by equation (2) may be applied in equalisation: 
the channel H(.) then represents the transmission channel 
in base band and X^: represents a series of symbols in the 
modulation alphabet, the width of which is equal to the 
channel memory. 
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During decoding of a convolutive code^ H(.) then 
represents the code generator polynomials^ and X;c 
represents the memory contents of the channel coder. 

We will now describe the concept of the likelihood 
5 ratio calculated in the context of the MAP algorithm in 
more detail. 

The likelihood ratio A^^ of a set X belonging to the 
set of parts of E is defined as being the ratio of the a 
posteriori probability of set X to the a posteriori 
10 probability of the complement of X in E. 

The likelihood ratio can be expressed in the form of 
the following equation: 

where PAP^ = l{Xj, =i /Yq ..yo...,..., Yk = Yk], 0 ^ k ^ K. 

15 Classically, the MAP algorithm is associated with a 

lattice representation of the different states to be 
considered. Thus, if a convolutive coder is considered, 
the MAP algorithm uses the lattice representation of the 
variation of the coder with time, for decoding. The 

20 lattice is usually represented by a graph, in which the 
abscissas axis represents the time and the ordinates axis 
represents the different possible states. Therefore each 
node in the lattice is associated firstly with a state, 
and secondly with an instant k. The possible transitions 

25 between nodes of two successive instants k and k+1 are 
represented by branches. 
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The MAP algorithm involves two intermediate 
variables, called the "forward variable" and the 
"backward variable" respectively. These variables 
propagated by the MAP algorithm are calculated 
5 recursively, in the forward direction and in the backward 
direction respectively with respect to the lattice, and 
the result is a calculation of the likelihood ratio at 
each instant k. 

These forward and backward variables defined in the 
10 remainder of the document use the probability of 
observing yk conditionally in the Xk = j and Xk-i = i 
states, expressed by the following equation: 

bi,j(Yi,) = Frob[Yi, = yn/Xi,.! ^ i, Xj, = j] (5) 

Classically, this probability is called metric. It 
15 may possibly depend on k, and is a function of the 
probability density of the noise B^. 

The MAP algorithm proceeds in three steps, firstly a 
forward recursion step, secondly a backward recursion 
step, and finally a data merging step. 
20 Forward recursion consists of calculating the 

probabilities af of passing through nodes representing 
the states Xi = i at instants k, starting from information 
obtained by passing through the lattice in the direct 

direction. The forward variable a* is thus defined as the 
25 probability of the event {Xi = i, Yo = Yor Yy^ = yjc) : 
af= Projb[Xi=i, Yo^^yor-^Yk^Yk] (6) 
It can be demonstrated that values of , 0 ^ k ^ K 
are calculated recursively as follows: 



with the initial value af = Vi, Vi G E. 

Similarly^ the backward recursion consists of 
calculating the probabilities /3f of passing through nodes 
representing the states Xi = i at instants starting 
from information obtained by going through the lattice in 

the indirect direction. The backward variable is thus 

defined as the probability of the event 
(Xi = i, Yfc-,1 = Yk-hi^ ... Yk = Yk) : 

i8,*= Prob[Xi = Yk^i = Yk^i, ... Yk = Yk] (7) 

It can be demonstrated that Pi , 0 ^ k ^ K are 
calculated recursively as follows: 

with the' initial value /?* = 1, Vi G E. 

The third step used by the MAP algorithm consists of 
merging the results of the two recursions to calculate 
the "a posteriori" probability of the states, by 
calculating the likelihood ratio of the detected symbols. 

The likelihood ratio of a set of states defined by 
equation (3) can then be expressed as a function of the 
forward and backward intermediate variables in the 
following form: 
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For example, in the context of symbol detection, the 
set X groups states corresponding to sending the same 
modulation symbol. In the context of channel decoding, 
the set X groups states corresponding to the same word at 
5 the input to the channel coder. 

Although the MAP algorithm supplies very high 
performance results, it has the disadvantage that it is 
complex to implement. 

In particular, the hardware implementation of the 
10 MAP algorithm is very difficult, and consequently 
receivers implementing this algorithm are expensive. 

Problems with implementation of the classical MAP 
algorithm are partly due to the strong dispersion of 

probabilities that therefore requires high precision 

15 for the calculation of the "forward" variable and the 

"backward" variable jSf . Moreover, calculation of the 

likelihood ratio as defined by ^equation (3) above means 

that divisions have to be imposed that are expensive and 

that can cause stability problems. 
20 Log-MAP and Max-Log-MAP type algorithms have been 

proposed to solve these various disadvantages with the 

MAP algorithm. 

Nevertheless, Log-MAP and Max-Log-MAP type 

algorithms, like MAP type algorithms, remain very complex 
25 to implement and require a large number of operations 

during the calculation of the likelihood ratio. 

In particular, one purpose of the invention is to 

overcome these disadvantages according to prior art. 
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More precisely, one purpose of the invention is to 
provide a signal processing technique implementing a MAP 
algorithm that is less complex than techniques according 
to prior art. More particularly, one purpose of the 
5 invention is to enable a reduction in the number of 
calculations used, both by MAP type algorithms, and by 
Log-MAP or Max-Log-MAP type algorithms, and achieving a 
better performance / complexity ratio than in prior art. 

Another purpose of the invention is to use such a 
10 technique, enabling the design of less complex receivers, 
and therefore less expensive than in the prior art. 

Another purpose of the invention is to provide such 
a technique that makes a satisfactory compromise between 
firstly low complexity requirements, and secondly 
15 reliability of the results. In other words, the invention 
is intended to provide a method of systematically 
reducing the complexity of MAP, Log-MAP and Max-Log-MAP 
type algorithms, with a good approximation quality. 

Thus, when classical MAP type algorithms usually 
20 necessitate the use of a number of operations in A^""^, 
where L represents the length of the transmission 
channel, and where A represents the number of characters 
in the modulation alphabet, the invention is intended to 
provide a technique for implementing a MAP type algorithm 
25 using a number of operations in a", where M < L-1, while 
maintaining a good approximation quality of the 
likelihood ratio. 

Another purpose of the invention is to implement 
such a technique that enables taking flexible decisions. 
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These objectives^ and others that will become clear 
later, are satisfied by implementing a signal processing 
method that uses a MAP (Maximum A Posteriori) type 

algorithm to determine a likelihood ratio of a set of 
5 states X of a lattice at instant each of the said 

states being associated with at least two intermediate 
variables called the "forward" variable and the 
"backward" variable respectively, propagated by the said 
MAP algorithm and calculated recursively, in a direct 

10 direction and in an indirect direction respectively at 
the said instant k with respect to the said lattice. 

According to the invention, such a method includes a 
step for reducing the number of states selected by the 
said MAP type algorithm during a calculation of the said 

15 likelihood ratio, and at least one determined value is 
assigned to the corresponding said forward and / or 
backward variable, so as to calculate an approximate 
likelihood ratio, at least for some non-selected states. 

Thus, the invention is based on a new and inventive 

20 approach to the use of a MAP algorithm that simplifies 
the different signal processing steps while maintaining a 
satisfactory quality of results. The invention proposes 
to reduce the number of states selected by the MAP 
algorithm to calculate the likelihood ratio, so as to 

25 simplify the different steps used during such a 
calculation, and to assign determined values to the 
missing intermediate variables in order to calculate an 
approximate likelihood ratio. Therefore, the invention is 
based on an approximate MAP type algorithm with a lower 
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calculation cost than conventional MAP algorithms. 
Therefore signal receivers using such an algorithm are 
less complex to construct than receivers according to 
prior art. 

5 Advantageously, for a given instant k, the said at 

least one determined value a(k) assigned to the said 

"forward" variable is such that 0 ^ a(k) ^ Min {af ) , and / 

or the said at least one determined value b(k) assigned 
to the said "backward" variable is such that 

10 0 ^ b(k) ^ Min { ) , where M/ and represent a set of 

the said states selected in the said direct direction and 
in the said indirect direction respectively at the said 

instant k, and where af and /3* represent the said 

forward and backward variables respectively at the said 
15 instant k. 

Preferably, for a given instant k, the said 
determined value a(k) and / or b(k) is unique and is 

assigned to at least one "forward" variable a* and / or 

"backward" variable . 

20 According to one advantageous embodiment of the 

invention, a constant value is assigned to the said 
"forward" and "backward" variables respectively, such 
that the said MAP type algorithm is a single-directional 
direct or indirect type algorithm respectively. 

25 According to one advantageous characteristic of the 

invention, the said step to reduce the number of states 
uses a "breadth-first" type lattice search algorithm. 
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According to a first preferred variant, the said 
"breadth-first" type algorithm is an M type algorithm. 

According to a second preferred embodiment, the said 
"breadth-first" type algorithm is a T type algorithm 
5 using at least one threshold. 

Preferably, the said at least one threshold is 
variable as a function of the said instant k. 

Advantageously, a predetermined value is assigned to 
the said variable threshold for each instant k. 
10 Preferably, for each instant k, the value of the 

said variable threshold is determined by the use of an 
adaptive algorithm. 

Advantageously, the said adaptive algorithm is a 
gradient type algorithm. 
15 According to one advantageous characteristic of the 

invention, when the said lattice comprises several nodes 
each associated with one of the said states and at a 
given instant k, the value of the said variable threshold 
T at an instant (k+1) is determined by the following 
20 equation: 

T(k-hl) = T(k) - u(M(k) - Mc) 
where T(k) represents the value of the said variable 
threshold at the said instant k, 

Mc is the target number of propagated nodes in the 
25 said lattice, 

lyi(k) is the number of propagated nodes in the 
said lattice at instant k, 

and 11 is a positive constant representing a 
learning gain. 



Preferably, the said adaptive algorithm is a 
gradient type algorithm with variable pitch. 

Advantageously, the said learning gain p is a 
function of the said instant k. 

Preferably, when the said "breadth-first" type 
algorithm is an M type algorithm, the said determined 
values a{k) and / or b(k) assigned to the said "forward" 
and / or "backward" variables respectively at a given 
instant k are given by the following equations: 

a(k) - Min (a^ )-Cf 
h(k) = Min (Pl')-ct 

where Cf and Ci, are two positive constants. 

Preferably, when the said "breadth-first" type 
algorithm is a T type algorithm, the said determined 
values a(k) and / or b(k) assigned to the said forward 
and / or backward variables at a given instant k 
respectively, are given by the following equations: 

a(k) = T^(k) - Cf 

h(k) = T^(k) - Ci, 

where Cf and are two positive constants, and where 
T^(k) and T^(k) denote the value of the said variable 
threshold at instant k in the said direct direction and 
in the said indirect direction respectively. 

Advantageously, the said MAP type algorithm belongs 
to the group comprising: 

- MAP type algorithms; 

- Log-MAP type algorithms; 

- Max-Log-MAP type algorithms. 



13 



It will be noted that the "forward" and "backward" 
variables may have different definitions depending on the 
type of algorithm used (MAP, Log-MAP or Max-Log-MAP) , as 
described in detail in the remainder of the document. 
5 Preferably, the said MAP type algorithm is a single- 

directional algorithm, and the said method uses a step to 
compare decisions made by the said single-directional 
algorithm with the corresponding decisions made by a 
Viterbi type algorithm, called Viterbi decisions. 

10 In this presentation, and in the remainder of the 

document, a Viterbi algorithm is an algorithm with a 
small number of states with which a "trace-back" 
operation is associated (in other words, decoding of the 
algorithm) . Thus, Viterbi decisions are decisions output 

15 by an algorithm with a reduced number of states, with 
which a "trace-back" is associated. 

According to a first advantageous characteristic, in 
the case of a negative comparison for at least one of the 
said decisions made by the said single-directional 

20 algorithm, the said method uses a substitution step for 
the said Viterbi decision corresponding to the said 
decision made by the said single-directional algorithm, 
called the substituted decision. 

Preferably, a determined value V is assigned to the 

25 absolute value of the said likelihood ratio associated 
with the said substituted decision. 

Advantageously, the said determined value V is equal 
to the absolute value of the average likelihood ratio of 
the sequence. 
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Therefore^ V corresponds to the average of 
reliabilities, which is approximately equal to the 
average reliability of the sequence. 

According to a second advantageous characteristic, 
5 in the case of a negative comparison for at least one of 
the said decisions made by the said single-directional 
algorithm, the said method uses a weighting step of the 
said likelihood ratio associated with the said decision 
considered, taking account of the said Viterbi decision. 
10 Preferably, when Y is a set of states associated 

with a decision output by the said Viterbi type 

algorithm at instant i, and represents the likelihood 
ratio associated with Y at instant i as calculated by the 
said single-directional algorithm during the said 

15 weighting step, the value of is replaced by A^ defined 

by A^ == A^ + X V, where V is a determined value. 

The invention is advantageously applicable to at 
least one of the domains belonging to the group 
comprising : 
20 - symbol detection; 

- signal coding / decoding; 

- turbo-decoding; 

- turbo-detection; 

- source coding by quantification (particularly 
25 quantification in lattice) . 

The invention also relates to a communication 
signals receiver comprising means of implementation of a 
MAP (Maximum A Posteriori) type algorithm to determine a 



likelihood ratio of a set of states X of a lattice at 
instant k, 

each of the said states being associated with at least 
one intermediate variable belonging to the group 
comprising a so-called "forward" variable and a so-called 
"backward" variable propagated by the said MAP algorithm 
and calculated recursively in a direct direction and in 
an indirect direction respectively at the said instant k 
with respect to the said lattice. 

According to the invention, such a receiver 
comprises means of reducing the number of states selected 
by the said MAP type algorithm in order to make a 
calculation of the said likelihood ratio, and for at 
least some non-selected states, at least one determined 
value is assigned to the corresponding said forward 
variable and / or backward variable, so as to calculate 
an approximate likelihood ratio. 

Other characteristics and advantages of the 
invention will become clearer after reading the following 
description of a preferred embodiment given as a simple 
illustrative and non-limitative example, and the attached 
drawings among which: 

Figure 1 shows a lattice representing a 
plurality of states associated with a MAP algorithm 
according to the invention; 

Figure 2 represents a flowchart illustrating 
the comparative or weighted comparative "trace-back" 
principle, used according to the invention in combination 
with single-directional MAP algorithms; 
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Figure 3 shows a block diagram of the different 
calculation steps used by the single-directional MAP type 
algorithm in Figure 2, applied to detection of symbols in 
an alphabet with A words, where A = 2^; 
5 - Figure 4 illustrates the different steps in the 

calculation of forward variables used by a two- 
directional MAP type algorithm according to the 
invention; 

Figure 5 .describes the different steps in the 
10 calculation of backward variables used by the algorithm 
in Figure 3, and the different steps in the calculation 
of likelihood ratios. 

The general principle of the invention is based on a 
reduction in the number of states selected during the 
15 calculation of the likelihood ratio of the MAP algorithm, 
so as to produce a simplified MAP algorithm. In one 
preferred embodiment of the invention^ this type of 
reduction is achieved particularly by using a "breadth- 
first" type algorithm. 
20 Before starting a detailed description of the 

figures mentioned above, we will present the principle of 
a reduction in the number of states used by the 
invention, and the different equations determined by the 
inventors . 

25 In one preferred embodiment of the invention, the 

inventors envisaged the use of "breadth-first" type 
algorithms, and particularly M type and T type algorithms 
to implement the MAP algorithms with a reduced number of 
states . 
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As for the Viterbi algorithm (described particularly 
in the article by G. D. FORNEY "The Viterbi algorithm". 
Proceedings of the IEEE, vol. 61, No. 3, March 1973), the 
M algorithm is a lattice "breadth-first" type search 
5 algorithm (this type of algorithm examines the cost 
associated with a family of states at every instant) . The 
M type algorithm is described particularly by J. 
B. Anderson and S. Mohan in "Sequential Coding 
algorithms; A survey and cost analysis", IEEE Trans, on 

10 Communications, Vol. COM-32, February 1984. 

For each lattice depth in Figure 1, the Viterbi 
algorithm only keeps the best path passing through each 
node, shown in thick lines in figure 1. 

The M type algorithm keeps the M best paths thus 

15 passing through the M best nodes, where M is an integer 
less than the number of nodes in the lattice of the 
Viterbi algorithm. Thus, when using the M algorithm, the 
nodes in the lattice (symbolised by black dots on the 
lattice in Figure 1) need to be classified as a function 

20 of their cumulated metric (note that the cumulated metric 
of nodes in the lattice is equivalent to the forward 
variable) so as to select the M best nodes. 

The T type algorithm is also a "breadth-first" type 
algorithm but is different from the M algorithm by the 

25 method used to select the nodes in the lattice in Figure 
1 propagated by the algorithm. The T algorithm selects 
only the nodes in the lattice associated with a cumulated 
metric for which the distance to the best cumulated 
metric of the lattice is less than a threshold T. 
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According to the invention, it is advantageously 
planned to introduce new "breadth- fir st " type T 
algorithms with a variable threshold. The inventors 
envisaged two distinct approaches to achieve this: 
5 - in a first variant embodiment of the invention, 

a threshold calculated a priori is assigned to each 
transition (remember that a transition from one state to 
another is symbolised by a branch of the lattice in 
Figure 1} . This threshold is calculated by using the 

10 finite length of the observation sequence (the 
observation sequence on the lattice in Figure 1 begins at 
instant k and terminates at instant k+4). Thus, if a 
lattice with length K+1 is considered, in which BMAX 
represents the maximum value of the branch metrics, the 

15 threshold T(k) will be chosen such that: 
T(k) ^ (K+l-k)BMAX and 0 ^ k ^ K; 

in a second variant embodiment of the 
invention, the value of the threshold of the T algorithm 
is fixed using an adaptive algorithm. An advantage of the 

20 T type algorithm compared with the M algorithm is that 
with a T type algorithm, there is no need to classify 
nodes in the lattice as a function of their cumulated 
metric in order to select the best nodes. The T type 
algorithm only imposes that the best node in the lattice 

25 is found in order to normalize cumulated metrics of all 
nodes in the lattice. In this second variant embodiment 
of the invention, the inventors have envisaged using an 
adaptive algorithm to fix the threshold of the T 
algorithm in order to slave its complexity. 
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For example, the threshold can be updated using a 
gradient algorithm and the following equation: 



where Mc represents the target number of propagated 
nodes, M(k) represents the number of nodes propagated at 
instant k, and p is a small positive constant that is 
called the learning gain. 

The algorithm can also be generalised to the case in 
which ]i depends on k, as is the case for gradient type 
algorithms with variable pitch. 

We will now present the implementation of these 
"breadth-first" type algorithms in more detail in the 
context of a reduction in the number of states of MAP 
type algorithms. 

The invention advantageously proposes to assign 
default values a (k) and b (k) to the "forward" and 
"backward" variables respectively at each instant k, so 
that MAP type algorithms can be used with fewer states. 

Let M{ and be the sets of states selected in 

the forward and backward directions respectively at 
instant k. 

If the state Xk = i forms part of the set M{ , then 
this state was selected at instant k. 

Since the "forward" variables are associated with 

probabilities, they must satisfy the normalization 



T(k-hl) = T(k) - i2(M(k)'Mc) 



(9) 



N 



condition 




It can be demonstrated that if this 
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condition is to be satisfied, then a (k) must be chosen 
such that: 

a (k) - — i — Vaf (10) 

It is quite obvious that the forward variables af 

5 for i ^ M( are not accessible (since the associated 

states are not selected, by definition of M( ) , and 
therefore it is only possible to determine the limits of 
the variation interval of a (k) . Assuming that the M 

states elements of M( have the M highest a posteriori 
10 probabilities (or APP) , a (k) is bounded by: 

0 ^ a(k) ^ Min (a^ ) (11) 

If a (k) is fixed to 0 regardless of k, it can be 
considered that the unpropagated states have a zero APP 
and do not participate in the likelihood calculation 
15 given by equation (8). In other words, this is equivalent 
to assigning an absolute confidence to the choice of the 
M propagated states. 

On the other hand, if a (k) is fixed at its maximum 
value, then minimum confidence is assigned to the choice 
20 of the M propagated states. 

The same reasoning can be applied to "backward" 
variables, so that b(k) is determined such that: 

0 ^ h(k) ^ Min ( fif ) (12) 



Equations (11) and (12) above may be refined in the 
25 context firstly of using an M type algorithm, and 
secondly a T type threshold algorithm. 
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When an M type algorithm is used, the classification 
of states selected in the direct (or "forward") direction 
and the indirect (or "backward") direction are used at 
each instant k, as a function of the values of af and 
5 respectively. Therefore, it is easy to determine the 
values of Mm(a^) and Mm^fil") , and therefore to assign the 

following values to a (k) and b(k): 

a(k) = Mm(a^) - Cf (13) 

b(k) = Mm(P^) - ct (14) 

10 where Cf and Cb are two positive constants. These two 
constants Cf and ct are used to adjust the confidence 
assigned to the choice of the propagated states. 

In the case of a T type algorithm, in other words 
when the states propagated by the algorithm are selected 

15 by a threshold (fixed or variable with time as proposed 
above) , the determination of the minimum values of the 
equations (11) and (12) requires a large additional 
calculation load. Consequently, for this class of 
algorithm, the inventors propose to fix a(k) and b(k) by 

20 means of the threshold: 

a(k) = T^(k)-Cf (15) 
b(k) = 'J^(k)-Cb (16) 
where T^(k) represents the threshold in the "forward" 
direction at instant k, and where (k) represents the 

25 threshold in the "backward" direction at instant k (for 
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example as defined by equation (9) above) , where Cb and Cf 
are two positive constants. 

In practice, Log-MAP or Max-Log-MAP type algorithms 
are frequently preferred to classical map algorithms, to 
5 simplify the implementation. The calculation of the 
likelihood ratio as defined by equation (3) above 
requires divisions for classical MAP type algorithms, and 
these divisions are expensive and may cause stability 
problems. This is why it is usually preferred to use 
10 algorithms working in the logarithmic range, that 
propagate the inverse logarithm of "forward" and 
"backward" variables and lead to the calculation of the 
logarithm of the likelihood ratio, also called Log- 
Likelihood Ratio (LLR) defined by: 



15 LLRf = -Log 



pap: 



(3a) 



where PAP'^ = P/Xjt = i / Yo = yo, Yk - YkJ^ 0 ^ k ^ K. 

If we consider a convolutive coder generating 
symbols X;c transmitted by a noisy channel to a receiver 
that receives symbols and implements a Log-MAP type 

20 algorithm to decode the received bits, the sign of LLRf 

gives the value of the decoded bit and its absolute value 

indicates the decision quality. 

Application of the reduction in the number of states 

principle presented above for classical MAP algorithms, 
25 to Log-MAP or Max-Log-MAP type algorithms is presented 

below: 
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Let a* = -Logia^) and fi] = -Log{pf) be the 
"forward" and "backward" logarithmic variables. According 
to the invention, the missing logarithmic variables for 
Log-MAP and Max-Log-MAP type logarithms with a reduced 
5 number of states are replaced in the forward and backward 

directions by a{k) and b{k) respectively, defined as 
follows : 

a{k) a: Mcixiai) (17) 

b{k) s Max(^]) (18) 

10 Thus according to the invention, when the reduction 

in the number of states of the Log-MAP or the Max-Log-MAP 
type algorithm is done using a "breadth-first" algorithm 
using a sort algorithm (in other words an M type 
algorithm), the missing data are replaced by: 

15 a{k) ^ Max{ai)+Cf for the forward direction (19) 

b{k) ^ Max{p))'\'Cb for the backward direction (20) 



where Cf and ct are two positive constants - 

Concerning "breadth-first" algorithms using a 

threshold to select states propagated in the "forward" 
20 and "backward" directions (in other words T type 
algorithms), the missing data are replaced by: 

a{k) = (k)+Cf for the forward direction (21) 

b{k) = (k)+Cb for the backward direction (22) 
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where T (k) and T (k) are the forward and backward 
thresholds respectively (possibly varying as a function 
of the time k, as proposed above) . 

The invention also provides a means of combining 
5 aspects related to simplification of classical MAP type 
algorithms (in other words MAP, Log-Map and Max-Log-MAP) 
with the principle of reducing states described above. 

As mentioned above, the MAP, Log-MAP and Max-Log-MAP 
algorithms require the propagation of forward and 
10 backward variables in the direct and indirect directions 
respectively, compared with the associated lattice. 

In one variant embodiment of the invention, it is 
proposed to simplify these algorithms by fixing "forward" 
and "backward" variables equal to a predetermined value, 
15 preferably constant, so as to produce single-directional 

k 

algorithms. Thus, the "forward" variables, a. and «/ can 

be fixed equal to an arbitrary constant value regardless 
of the values of k and i, to obtain "indirect" MAP, Log- 
MAP or Max-Log-MAP algorithms. Similarly, the backward 

20 variables or can be fixed equal to an arbitrary 

constant value regardless of the values of k and i, to 
obtain "direct" algorithms. 

The principle of reducing the number of states can 
then be applied to these direct and / or indirect 

25 algorithms using "breadth-first" type algorithms 
described above, so as to obtain direct and / or indirect 
algorithms with a reduced number of states. Thus, for 
example, in the case of an indirect MAP algorithm for 
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which all forward variables a* have been fixed to a 

constant value, the missing backward variables 
(corresponding to unselected states) are replaced by a 
predetermined value h(k), as described above in the 
5 document . 

This type of missing data management makes direct 
and indirect algorithms using a small number of states 
more efficient. When the number of propagated states is 
small, it is probable that at least one of the sets on 

10 which the likelihood ratio defined by equation (3) is 
calculated, will be empty. 

The proposed method, according to which default 
values are assigned to the intermediate variables 
necessary for the calculation of the likelihood ratio but 

15 corresponding to unpropagated states of the lattice, can 
then be used to determine the likelihood ratio as soon as 
there are at least two states propagated in the lattice, 
and to assign a soft value to it. If the value 0 is 
assigned to the probability of an unpropagated state, the 

0 

20 values of the likelihood ratio equal to or 

constant 

constant 

will be obtained for MAP type algorithms. For 

0 

Log-MAP type algorithms, values of the likelihood ratio 

equal to -oo or -hoc are obtained. 

One disadvantage of these single-directional 
25 ■ algorithms that has been identified by the inventors, is 
that decisions made by these algorithms^ in other words 
the sign of the likelihood ratios defined by equation (3) 
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(or (3a) ) are usually different from decisions taken by 
algorithms using a " trace-back" . (Remember that 
algorithms using a "trace-back" enable a decision to be 
taken using the most likely path) . 
5 The transition matrix Yl flight not allow the states 

associated with decisions made by the different single- 
directional algorithms. Therefore, in one new 
advantageous embodiment of the invention, it is envisaged 
to use a single-directional algorithm with a small number 

10 of states with a comparative trace-back or a weighted 
comparative trace-back. 

According to the invention, the comparative "trace- 
back" consists of comparing decisions made by a single- 
directional algorithm with the corresponding decisions 

15 made by a Viterbi type algorithm as shown in Figure 2 . A 
set of states Y is considered, associated with the 

decision Df output by the Viterbi type algorithm at time 

i. Let Af be the likelihood associated with Y at time i, 
calculated by the single-directional algorithm, and let 
20 Af be the corresponding decision taken by the single- 
directional algorithm. 

During a first step reference 21, the likelihood 

ratio Af is calculated so as to make a decision Af for 
the single-directional algorithm. The corresponding 
25 decision of for the Viterbi type algorithm is analysed 
(22) and these two decisions Af and of are compared (23) . 
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If the two decisions are different, then the 
decision taken by the Viterbi type algorithm replaces 
(25) the decision made by the single-directional 
algorithm. The absolute value of the likelihood ratios 
5 assigned to these decisions remains unchanged or is 
forced to a value V. V is fixed a priori to be equal to a 
small value, or is calculated for each sequence. If it is 
calculated for each sequence, V may for example be fixed 
equal to the absolute value of the average likelihood 
10 ratio of the sequence. 

Rather than having to make a clear choice between 

the decision made by the Viterbi type algorithm Df and 
the decision produced by the single-directional algorithm 
Af, it would also be possible to weight these two 
15 decisions: this is called the weighted comparative "trace 
back" principle, also shown in figure 2. 

When the comparative weighted "trace back" step has 

been completed, the likelihood Af is replaced (24) by 

A^ defined by: 

20 A^ = Af + Df X V (23) 

where V is equal to the absolute value of the average 
likelihood ratio calculated on the sequence, or fixed a 
priori . 

We will now present an embodiment of the different 
25 calculation steps used by a single-directional 
logarithmic iVIAP type algorithm according to the 
invention, with reference to Figure 3, applied to 
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detection of symbols in an alphabet with A words, where A 
= 2^ 

The CM(i) vector 31 contains cumulated metrics M(i) 
corresponding to the M(i) states retained at instant i 
5 stored in the State (i) vector 32 with size M(i) times Q 
bits . 

The "calculate BM and ACS" block 33 calculates 
A.M(i) cumulated metrics associated with the A.M(i) 
states following the M(i) states in the State (i) vector 

10 32. It then selects the best states converging towards 
the same node. This selection consists of either 
selecting the states that have the smallest cumulated 
metric (in the case of a single-directional Max-Log-MAP 
type algorithm) ^ or recursively performing a correction 

15 comparison operation by using the Jacobian logarithm 
formula (in the case of a single-directional Log-Map type 
algorithm) . 

Thus, in the case of a single-directional Log-Map 
type algorithm, this method includes three main steps: 
20 a) Step 1: calculate the branch metric BM(i^j) = 

-ln(bij(yk)) in accordance with equation (5) described 
above. Note that the branch metric is expressed simply as 
a Euclidean distance. 

b) Step 2: accumulation. 

25 The variables (respectively /8f ) and are updated 

in a simplified manner. At instant k-1, two nodes Po and 
Pi are considered associated with the forward variables 

^m/) ^^u/) respectively. At the next instant k, the 
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final node associated with the cumulated metric Ciyij(k) is 
considered. The branch metric from node Po to the final 
node is BM(Po(j)f j) and the branch metric from node Pi to 
the final node is BM(Pi(j), j). We then get: 

5 a* = maK[{a'p-l^-hBM(Pi(j) , (a'^\jy'tBM(Po(j),j))] 

c) Step 3: selection 

The best path that converges towards each node is 
selected, and this step is used for the "trace back". 

The first and third steps for a Log-MAP type 
10 algorithm are identical, but during the second step to 
update the forward variables, we use: 

aJ = max[{a'jr,lj^+BM(Pi(j) ,j) ) , (cl'-^^j^+BM (Pq (j ) r j ) ) ] + 

f(\{CL'-^U)^BM(Pi(j),j) - (o,';^^j^-hBM(Po(j) rj) )\, 
where f(z) = ln(l + exp{-z)). 

15 The metrics of the best nodes selected by the 

"calculate BM and ACS" block 33 are saved in the 
TmpCM(i+l) vector 34. The number of states selected in 
this step is between M(i) and A.M(i). The corresponding 
states selected in this step are stored in the 

20 TmpState (i+1) vector 35. 

The "select M or T" block 36 is used to select 
states propagated by the MAP algorithm. 

If the reduction in the number of states uses a T 
type algorithm, the "select M or T " block 36 selects 

25 states for which the cumulated metric distance starting 
from the best cumulated metric is less than a threshold 
T. 
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If the reduction in the number of states uses an M 
type algorithm, the "select M or T" block 36 selects the 
• best M states. 

In the case of a T type algorithm, the "select M or 
5 T" block 36 also calculates the new threshold T(i+1) as a 
function of the number of selected states M(i), applying 
the threshold T(i) as described previously in equation 
(9) . 

After the selection step performed by the block 

10 reference 36, the M(i+1) cumulated metrics corresponding 
to the M(i+1) states selected at instant i+1 stored in 
the State (i+1) vector 38 with size M(i+1) times Q bits, 
are stored in the CM{i+l) vector 37. 

The Survivor { i+l ) vector 39 with size equal to N Q- 

15 bit words contains the index of the M(i) branches that 
converge towards the selected M(i) states stored in the 
State (i) vector 32. The state towards which the winning 
branch of the lattice converges is coded by the address 
at which the branch index is stored, on Q bits. 

20 The "trace back" operation is performed using all 

Survivor (i+l) vectors 39, for i varying from 0 to K-1, in 
order to determine decisions to be compared with or 
weighted by decisions determined by the LLR (Log- 
Likelihood Ratio) calculation block 40, as described 

25 above in relation to Figure 2. 

We will now describe an embodiment of the forward 
variables calculation for a two-directional MAP logarithm 
according to the invention, in relation to Figure 4. 



r 
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Remember that two-directional Log-MAP or Max-Log-MAP 
algorithms with a small number of states begin with a 
first step in which the forward variables are calculated 
recursively in the same way as single-directional 
5 algorithms. 

Figure 4 illustrates the calculation of the list of 
states kept and cumulated metrics in the direct or 
forward direction. 

Unlike single-directional algorithms presented in 
10 Figure 3, the cumulated metrics of the CMF(i} vector 42, 
the corresponding selected states of the State F(i) 
vector 43, and possibly the number of states selected at 
each instant, are stored so that they can be used to 
calculate the likelihood ratios LLR during the backward 
15 step. 

As shown in Figure 3, the CMF(i) vector 42 contains 
the cumulated metrics MF(i) corresponding to the iyiF{i) 
states selected in the forward direction at instant i 
stored in the StateF(i) vector 43 with size equal to 
20 MF(i) times Q bits. 

The "calculate BM and ACS" block 33 calculates the 
A.MF(i) cumulated metrics associated with the A.MF(i) 
states following the MF(i) states of the StateF(i) vector 
43. Then, it selects the best states converging towards 
25 the same node, using a method similar to that described 
with reference to Figure 3. 

The metrics of the best nodes selected by the 
"calculate BM and ACS" block 33 are stored in the 
TmpCMF(i+l) vector 44. The number of states selected in 
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this step is between MF(i) and A.MF(i). The corresponding 
states selected in this step are stored in the 
TmpStateF (i+1) vector 45. 

The "select M or T" block 36 is used to select 
5 states propagated by the MAP algorithm. 

If the reduction in the number of states uses a T 
type algorithm, the "select M or T" block 36 selects the 
states for which the distance between the cumulated 
metric and the best cumulated metric is less than a 
10 threshold T. 

If the reduction in the number of states uses an M 
type algorithm, the "select M or T" block 36 selects the 
M best states. 

In the case of the T type algorithm, the "select M 
15 or T" block 36 also calculates the new threshold TF(i+l) 
in the forward direction as a function of the number of 
selected states iyiF{i) by applying the threshold TF(i) as 
previously described in equation (9). 

After the selection step carried out by the block 
20 reference 36, the MF(i+l) cumulated metrics corresponding 
to iyiF(i+l) states selected at instant i + 1 stored in the 
StateF(i+l) vector 47 with a size equal to MF(i+l) times 
Q bits, are stored in the CiyiF(i+l) vector 46. 

The cumulated backward metrics of the CMB(i) vector 
25 50 and the selected states stored in the StateB(i) vector 
51 are calculated symmetrically to the forward 
calculation as shown in Figure 5. Figure 5 also presents 
the calculation of LLR values from forward variables 
(cumulated metrics, list of states and possibly 
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thresholds) memorised in the direct direction, and new 
variables calculated in the indirect direction. 

Thus, the CMB(i) vector 50 contains the MB(i) 
cumulated metrics corresponding to the MB(i) selected 
5 states in the backward direction at instant i stored in 
the StateB(i) vector 51 with size equal to Mb(i) times Q 
bits . 

The "calculate BM and ACS" block 33 calculates the 
A.MB(i) cumulated metrics associated with the A.MB(i) 
10 states following the MB(i) states of the StateB(i) vector 
51. The next step is selection of the best states 
converging towards the same node, using a method similar 
' to that described with reference to Figures 3 and 4. 

The metrics for the best nodes selected by the 
15 "calculate BM and ACS" block 33 are saved in the 
TmpCMB ( i+1 ) vector 53. The number of states selected in 
this step is between MB(i) and A.MB(i). The corresponding 
states selected in this step are stored in the 
TmpStateB (i+1) vector 52. 
20 The "select M or T" block 36 is used to select the 

states propagated by the MAP algorithm. 

If the reduction in the number of states uses a T 
type algorithm, the "select M or T" block 36 selects the 
states for which the distance between the backward 
25 cumulated metric and the best cumulated metric is less 
than a threshold T. 

If the reduction in the number of states uses an M 
type algorithm, the "select M or T" block 36 selects the 
M best states. 
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In the case of a T type algorithm, the "select M or 
T" block 36 also calculates the new. threshold in the 
backward direction Tb(i+1) as a function of the number of 
selected states MB(i) by applying the TB(i) threshold as 
5 described above in equation (9) • 

After the selection step performed by the block 
reference 36, the MB(i+l) cumulated metrics corresponding 
to the MB(i+l) states selected at instant i+1 stored in 
the StateB(i + l) vector 55 with size equal to MB (i + 1) 
10 times Q bits, are stored in the CMB(i+l) vector 54. 

The LLR calculation 40 uses the cumulated metrics of 
the CMF(i) vector 42 and the CMB(i) vector 50, and the 
selected states stored in the StateF(i) vector 43 and the 
States (i) vector 51. For the forward direction, the 
15 cumulated metrics of the CMF(i) vector 42, the states of' 
the StateF(i) vector 43, and possibly the thresholds 
TF{i) and the number of selected states MF(i) are stored 
for all sequences, in other words for i varying from 1 to 
K. 

20 The "calculate LLR" block 40 may use TF(i) and TB(i) 

thresholds in the case of a reduction in the number of 
states using a T type algorithm, or the less good 
cumulated metrics CMB(M) and CMF(M) in the case of a 
reduction in the number of states using an M type 

25 algorithm, to manage the missing data, as described 
above. In the same way as for single-directional 
algorithms, the LLR calculation uses the Max-Log-Map 
approximation or the Jacobian logarithm formula for Log- 
Map type algorithms. 



